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ABSTRACT 

We present photometric and spectroscopic observations of a peculiar variable (designated DASCH 
J075731. 1+201735 or J0757) discovered from our DASCH project using the digitized Harvard College 
Observatory archival photographic plates. It brightened by about 1.5 magnitudes in B within a year 
starting in 1942, and then slowly faded back to its pre-outburst brightness from 1943 to the 1950s. 
The mean brightness level was stable before and after the outburst, and ellipsoidal variations with a 
period of P = 119.18±0.07 days are seen, suggesting that the star is tidally distorted. Radial-velocity 
measurements indicate that the orbit is nearly circular (e = 0.02 ± 0.01) with a spectroscopic period 
that is the same as the photometric period. The binary consists of a 1.1 ± 0.3 Mq MOIH star, and 
a 0.6 ± 0.2 Mq companion, very likely a white dwarf (WD). Unlike other symbiotic binaries, there is 
no sign of emission lines or a stellar wind in the spectra. With an outburst timescale of ~ 10 years 
and estimated B band peak luminosity Mb ~ 0.7, J0757 is different from any other known classic or 
symbiotic novae. The most probable explanation of the outburst is Hydrogen shell-burning on the 
WD, although an accretion-powered flare cannot be ruled out. 
Subject headings: binaries: symbiotic — novae, cataclysmic variables 



1. INTRODUCTION 

Symbiotic stars are interacting binaries consisting of 
an evolved red giant and a hot object, which is usually 
a White Dwarf (WD; see e.g. Kenyon 1986, 1994; Nuss- 
baumer 2000; Mikolajewska 2007 and references therein). 
Mass is transferred from the red giant to its hot com- 
panion via stellar wind in most cases, but could also be 
Roche lobe overflow in some cases. They usually show 
low temperature absorption features from the red giant, 
and strong emission lines from surrounding circumstellar 
material ionized by the hot component. Most symbiotic 
stars have red giant masses 0.6 — 3.2 Mq, hot component 
masses 0.4 - 0.8 M©, and orbital periods ^ 200 - 1000 
days (Mikolajewska 2003). 

Symbiotic stars are highly variable sources with vari- 
ation timescales from minutes to decades. The most 
common type of outburst, called a classical symbiotic 
outburst, recurs multiple times on decades-long time 
scales with amplitudes of 1 — 3 mag, and with each out- 
burst lasts from months to a few years (Kenyon 1986; 
Mikolajewska 2011). The nature of these outbursts is 
not entirely clear, and may be related to instabilities 
in the accretion disks (Mikolajewska 2003; Sokoloski et 
al. 2006). Another relatively rare type of outburst is 
symbiotic nova (also called slow nova), which shows a 
single long outburst lasting from years to decades due 
to thermonuclear runaways at the surface of the WD 
(Allen 1980; Paczyhski & Rudak 1980; Kenyon & Web- 
bink 1984). The peak luminosities range from 4000 to 
40,000 Lq, depending on the mass of the WD (Miirset 
& Nussbaumer 1994; Iben 2003). They mimic classical 
novae but with much longer timescales. Only 9 sym- 
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biotic novae are known so far (Belczynski et al. 2000; 
Mikolajewska 2011). 

DASCH (Digital Access to a Sky Century @ Harvard) 
is a project to digitize and analyze the scientific data 
contained in the Harvard College Observatory (HCO) 
photographic plates taken from the 1880s through the 
1980s. The motivation is to explore the temporal varia- 
tions of stars and active galaxies on the relatively poorly 
explored ^ 1 — 100 yr timescales (but extending down to 
~ 1 week timescales as well), as summarized by Grind- 
lay et al. (2009 and 2012, in preparation). We developed 
the astrometry and photometry pipeline (Laycock et al. 
2010; Los et al. 2011; Servillat et al. 2011; Tang et al. 
in prep.), and have scanned ~ 19,000 plates in five ini- 
tial fields. Here we report the discovery of an unusual 
10-yr outburst in a peculiar symbiotic system, DASCH 
J075731.1-F201735 (named after its equatorial J2000 co- 
ordinate; hereafter J0757). 

2. OBSERVATIONS AND RESULTS 

2.1. DASCH and ASAS light curves 

J0757 was measured on 694 plates near the M44 field 
noted by its peculiar long-term variability. These plates 
cover 5—25 degrees on a side with typical limiting mag- 
nitudes B ^ 14—15 mag, and most of them are blue sen- 
sitive emulsions. We used the GSC2.3.2 catalog (Lasker 
et al. 1990) for photometric calibration, and our typical 
uncertainty is ^ 0.1 — 0.15 mag (Laycock et al. 2010; 
Tang et al. in preparation). There are ~ 1.2 x 10^ ob- 
jects with more than 100 magnitude measurements dis- 
tributed over ^ 100 years from these plates covering the 
M44 field. J0757, which is N2211021132 in the CSC cat- 
alog, is the only one found with > 1 magnitude outburst 
on a ~ 10 year timescale above a well-defined quiescence 
level. 

The DASCH light curve of J0757 is shown in Figure 
1. We supplement the figure with V-band data from 
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ASAS starting in the year 2000 (J0757 has the desig- 
nation ASAS J075731+2017.6; Pojmanski 2002). J0757 
is classified by ASAS as a semi-detached/contact binary 
with best-fit period of 119.2 days with 0.16 mag varia- 
tions in the V band. Since ASAS data are in the V band, 
while DASCH magnitudes are in B, we added 1.5 mag 
to the ASAS V mag in the plot, which is about the sum 
of the maximum galactic extinction {E{B — V) = 0.06; 
Schlegel et al. 1998) and the typical B — V value for 
a MOIII star (1.43; Pickles 1998). We will discuss the 
spectral classification in section 2.2. 

As shown in Figure 1, J0757 was in quiescence until 
about March 1942. J0757 then brightened over the rest 
of the year of 1942, and slowly decayed until about 1950. 
The peak magnitude was reached around January, 1943. 
There is a seasonal gap between April and October in 
1942 when the star was low or not visible in the night 
sky, and thus we do not know when exactly the out- 
burst started. If we naively assume that its magnitude 
increased linearly, then the outburst rise time is ^ 88.5 
days (see the lower right panel in Figure 1). 

Before and after the outburst when data are available 
(i.e. 1890-1942 and 1962 - 1990), the magnitudes of 
J0757 were reasonably stable (note that the relatively 
constant B~ 11.9 from ^ 1950 — 1954 was not yet the 
quiescent value) . The median DASCH magnitude before 
the outburst (e.g. before June 1942) is 12.24 and the me- 
dian magnitude after the outburst (after 1960) is 12.26. 
The median magnitude of the entire DASCH light curve 
excluding the period between June, 1942 and the start 
of 1960 is 12.25, with rms scatter = 0.14 mag for the 607 
data points outside the outburst. The median ASAS V 
magnitude is 10.75 mag, which corresponds to B=12.24, 
assuming the typical colors of an MOIII star with galactic 
extinction. 

2.2. Light curve folding and the photometric period 

Both the DASCH and ASAS light curves show ellip- 
soidal variations. The ellipsoidal variation amplitude is 
- 0.11 mag in the DASCH data (B band, 1890-1942 
and 1962 - 1990), and - 0.16 mag in the ASAS data (V 
band, 2003 - 2009). 

We folded the light curves of J0757 to search for its 
photometric period. We divided DASCH light curve into 
two phases, i.e. before the outburst and after the out- 
burst, as described in section 2.1. We then got four sets of 
light curve data: ASAS light curve, DASCH light curve 
excluding the outburst (before or after the outburst), 
DASCH light curve before the outburst, and DASCH 
light curve after the outburst. 

We first folded the ASAS light curve with a given pe- 
riod, and fit it with a 3 harmonics Fourier model (an 
example is shown as the blue line in the upper panel of 
Figure 2). We fixed the model line in shape, but allow 
two free parameters, i.e. median magnitude and ampli- 
tude, for adjustments. Next, we fit the four sets of light 
curve data to the 3 harmonics Fourier model, and derived 
the reduced as a function of different trial periods, as 
shown in the lower panels in Figure 2. DASCH data were 
binned by phase at the given period in the fitting. Our 
best fit period is 119.18 ± 0.07 days. The fact that the 
binary period is unchanged (within the errors) sets limits 
on the change of the orbit (due to mass transfer or onset 
of common envelope evolution). The resulting value for 



a/a is 0±4 X 10"^ yr^^. 

2.3. Optical Spectroscopy 

Spectra were obtained with the FAST spectrograph 
(Fabricant et al. 1998) on the 1.5-m Tillinghast refiec- 
tor telescope at the F. L. Whipple Observatory (FLWO) 
for spectral classification. The spectra were reduced and 
wavelength-calibrated with standard packages (Tokarz & 
Roll 1997). Figure 3 shows the FAST 300 grating spec- 
trum, which has a resolution of about 7 A. The upper 
panel shows the whole spectrum, which suggests that its 
spectral type is MO; The lower panels show two luminos- 
ity indicators, which suggest that DASCH J0757 belongs 
to luminosity class HI. The lower-left panel shows Ca I 
at 4226 A, which disappears with increasing luminosity, 
indicating that DASCH J0757 is not a luminosity class 
I star. The lower-right panel shows the CaH B-band 
at 6385 A, and the blend near 6362 A consisting of Ti, 
Fe, and Cr. The weakness of CaH B-band suggests that 
DASCH J0757 is not a main sequence star. We there- 
fore adopt a spectral type of -MOIII for DASCH J0757, 
with uncertainties of ±1 in both spectral type and 
luminosity class. 

We also obtained spectra using the MIKE echelle spec- 
trograph (Bernstein et al. 2003) on the 6.5 m Maga- 
Uen/Clay telescope with a 0.7" x5" slit on January 17, 
2011. The spectral resolution is i? ~ 30,000. The spec- 
tra were extracted using an IDL pipeline developed by 
S. Buries, R. Bernstein, and J. S. ProchaskgQ- The re- 
duced spectra have S/N ~ 50 - 100 in 4300-9000 A 
and S/N >~ 10 in 3700-4300 A. We computed a grid 
of spectrum templates for a range of abundances, effec- 
tive temperatures, and gravities using Kurucz's programs 
ATLAS12 and SYNTHE (Kurucz 2005). We derived the 
projected rotation velocity of 10 ± 1 km s^^ . In M giants 
the continuum opacity is mainly H- which varies with 
metal abundance. The ratio of line to continuum does 
not change strongly with abundance. The TiO back- 
ground is strong at any abundance because both Ti and 
O are alpha-enhanced at low abundances. We estimated 
an effective temperature of Te// = 3850 ± 50 K; a sur- 
face gravity of log g = 1.0 ± 0.5 cgs; and a metallicity of 
[Fe/H] = —0.6 ± 0.1, assuming J0757 is a population II 
star with alpha enhancements. The spectra can also be 
fit with Population I solar abundance templates. There 
is no sign of emission lines or wind in its spectra. 

To measure the radial velocity curve of J0757, we 
obtained spectra on 15 different nights from October, 
2010 through April, 2011 with the Tillinghast Reflec- 
tor Echelle Spectrograph (TRES; Szentgyorgyi & Furesz 
2007) on FLWO 1.5 m. These spectra were reduced 
and wavelength-calibrated with standard packages (Mink 
2011). The spectral resolution is i? ~ 44,000. The re- 
duced spectra have typical S/N ~ 10 - 30 in 4500-9000 
A, and S/N - 5 - 10 in 4000-4500 A. 

Radial velocities of DASCH J0757 were derived from 
the TRES spectra, and were fitted to a model velocity 
curve using the derived photometric period (see Figure 
4). The derived orbital parameters are listed in Table 1. 

http: //web.mit .edu/~burles/www/MIKE/ 
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The resulting mass function is 

/i(M2) = (3.31 ± 0.02) X IQ-^ Mq, (1) 

and 

Ma sin i = (0.3211 ± 0.0006) (Mi + Ma)^/^ M©, (2) 

where subscript 1 denotes the red giant, and 2 denotes 
its companion. 

For comparison, we have also fit a model velocity curve 
without using the photometric period of P = 119.18 ± 
0.07 days, and got similar results: Tq = 2455597.969 ± 
0.06, P = 118.89 ± 0.080 days, 7 = -33.57 ± 0.04 km 
s-i, K = 13.885 ± 0.023 km s-\ e = 0.024 ± 0.015, 
and w = 70.38 ± 0.09 deg. In the following discussion 
throughout the paper, we adopt the parameters obtained 
using the photometric period of P = 119.18 ± 0.07 days, 
as listed in Table 1. 

2.4. SED from archival data 

Using archival data from WISE, IRAS, AKARI/IRAC, 
2MASS, FON, SDSS, GALEX and ROSAT, we plot- 
ted the spectral energy distribution (SED) of J0757, as 
shown in Figure 5. The phase of the GALEX observation 
is 0.603 ± 0.01, and the phase of ROSAT observation is 
0.426 ± 0.02. Here phase 0.5 is defined as inferior con- 
junction of the red giant. The SED is well matched by 
a synthetic spectrum of a Tcff = 3750 K, and \ogg = 1.5 
giant (Lejeune et al. 1997), from IR to optical, as shown 
in black line in Figure 6. However, there is an excess in 
the GALEX NUV (~ 1900 - 2700A), which suggests a 
hot component in addition to the giant. 

3. SYSTEM PARAMETERS 
3.1. Ellipsoidal models 

We derive system parameters from the light and veloc- 
ity curves. For this purpose we used the ELC code of 
Orosz & Hauschildt (2000) to model the ASAS V-band 
light curve and the TRES radial velocity curve. The 
ELC model is based on standard Roche geometry and 
uses specific intensities from model atmospheres, in par- 
ticular the NextGen models of Hauschildt et al. (1999). 
ELC's use of the model atmosphere specific intensities 
eliminates the need for a parameterized limb darkening 
law, which is important in cases with red giants where 
the limb darkening near the stellar limb is very nonlinear 
(Orosz & Hauschildt 2000). We don't have model atmo- 
sphere specific intensities appropriate for white dwarfs. 
However, since the red giant dominates the V-band flux, 
this is not a serious limitation. 

In our initial runs, we had 8 free parameters: the or- 
bital period P, the time of periastron passage Tq, the 
inclination i, the Roche lobe filling factor /i, the eccen- 
tricity e, the argument of periastron a;, the i^T-velocity of 
the red giant and the mass of the red giant Mi. We 
assumed synchronous rotation. The models are some- 
what insensitive to the temperature of the red giant, and 
that parameter was held fixed at Ti = 3720 K. In addi- 
tion to the light and velocity curve, we used the projected 
rotational velocity of the red giant as an extra constraint. 

The model fits were optimized using ELC's genetic 
code and its Monte Carlo Markov Chain code. After in- 
spection of the post-fit residuals, additional free parame- 
ters were considered. We allowed the red giant's rotation 



to be asynchronous with the orbit, where 51 is the ratio 
of the spin frequency to the orbital frequency. We also 
accounted for heating by the WD. In this case, we fixed 
the luminosity of the irradiating source to logi^ = 36.6 
(c.g.s. units) and allowed the albedo Ai to be a free pa- 
rameter (the choice of Ai is somewhat arbitrary since dif- 
ferent combinations of Li and Ai can give identical light 
curves). Finally, we included an accretion disk that is 
parameterized by the outer radius in terms of the Roche 
radius rout, the temperature at the inner rim T^, and the 
exponent of the disk temperature distribution ^, where 
r(r) =Td(r/r,)«. 

Some typical results are shown in Figure 6. The binned 
ASAS light curve has maxima with different heights, 
which could be indicative of an eccentric orbit. How- 
ever, the simultaneous fit to the radial velocities indi- 
cates a small eccentricity (0.012). Also, the minimum 
near phase 0.5, corresponding to the inferior conjunction 
of the red giant, is deeper than the minimum at phase 
0.0. For a pure ellipsoidal light curve, the minimum at 
phase 0.0 would be the deeper one. Heating by the ir- 
radiating source "fills in" the minimum near phase 0. 
However, in all of the models considered, we never had a 
case where the model curve had a deeper minimum near 
phase 0.5, so simple heating cannot explain the difference 
in the observed minima. 

An accretion disk could potentially eclipse the red gi- 
ant near phase 0, thereby making the minimum there 
deeper. This is opposite of what we want. Normally, the 
accretion disk does not give phase modulated light in the 
ELC model (if it is not eclipsed by the star). However, 
a "hot spot" on the disk rim can be used to mimic the 
place where the matter transfer stream impacts the disk. 
As the binary turns in space, different areas of this spot 
are visible, giving rise to extra light at certain phases. 
An example model is shown in Figure 6. The hot spot 
has an angular width of 18°, an azimuth (measured from 
the line of centers in the direction or the orbital motion) 
of 66°, and has a temperature of 3800 K, which is 4.3 
times hotter than the rest of the disk rim. The maxi- 
mum near phase 0.75 is fit, but the minimum near phase 

0. 5 is still not well-fit. 

We did two runs of fitting with a Monte Carlo Markov 
chain, one with an external constraint on the rotational 
velocity, and one without the constraint. Results from 
both runs are similar. The best-fit ellipsoidal models 
have incfination of 60° ± 5°, WD mass of 0.5 ± 0.1 Mq, 
red giant mass of 1 ± 0.3 Mq, and red giant radius of 
43±4P0. 

We note that in order to match the observed rota- 
tional velocity of Vrotsini = 10 ± 1 km s~^ derived 
from the MIKE spectra, the red giant must be rotat- 
ing 0.7 times slower than its synchronous velocity. Given 
the typical tidal synchronization timescale of ~ 10^ yr 
(Zahn 1989) adopting the parameters listed in Table 

1, the red giant is expected to be co- rotating with the 
binary. Consequently, the size of the giant should be 
Pi = PFrot/27r = 23.5/ sin i Rq, which is 27 P© if 
i = 60°, which are smaller than the radius derived from 
ellipsoidal fitting. Such a discrepancy has been seen 
in many symbiotic binaries, and is an unresolved "em- 
barrassing problem" (Mikolajewska 2007). It probably 
comes from the extended atmosphere, and that limb 
darkening near the limb is very nonlinear. The light 
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curve models with near-linear limb darkening produce el- 
lipsoidal light curves with amplitudes that are too small. 
Also, the rotational broadening kernel can be a lot differ- 
ent than the usual analytic kernel, which in turn leads to 
strange biases in the Vrotsini value (Orosz & Hauschildt 
2000). 

3.2. Estimates from the spectroscopic observations 

Here we estimate the astrophysical parameters of 
DASCH J0757, based on the results in Sections 2 and 3.1. 
Since the galactic extinction is small {E{B — V)= 0.06; 
Schlegel et al. 1998), we do not take it into account for 
simplicity. The following estimated parameters are sum- 
marized in Table 1. 

The inclination from the ellipsoidal fitting is ?' ^ 55° — 
65°, and we adopt a more tolerant value of i = 45° — 75° 
in the following estimates. The radius of the red giant 
ranges from ~ 25 — 33 Rq (from Vrot assuming synchro- 
nized rotation) to ^ 43 Rq (from ellipsoidal fitting), 
as discussed in Section 3.1. Therefore, we adopt a ra- 
dius and error bar that encompasses both methods as 
i?i = 35 ± 9 i?0. At Teff = 3850 ± 50 K, the bolometric 
luminosity of the red giant is L = 250 ± 130 Lq, and its 
distance is 1.0 ± 0.3 kpc. The absolute magnitude of the 
system in quiescence is then Mb = 2.2 ± 0.6. 

Combining the radius, effective temperature, log g = 
1.0±0.5 cgs and [Fe/H] = — 0.6±0.1, the red giant mass 
inferred from stellar evolutionary models (Bertelli et al. 
2008) is Ml ~ 0.9 ± 0.4 Mq. This mass estimate for the 
giant is determined solely from the equations of stellar 
structure and is therefore independent of its prior evolu- 
tion. The mass of the giant may have been greater in the 
past. However, considering the giant is only partially fill- 
ing its Roche lobe as inferred from the ellipsoidal fitting, 
and given the MOIII wind loss rate of ~ 10~^ - 10~^ Mq 
yi~^ (Cox 2000), we do not expect it has lost a significant 
fraction of its current mass over the ~ 10^ yr lifetime of 
the giant in its current phase of evolution. Requiring 0.8 
Mq to turn-off from main sequence within a Hubble time, 
we then have Mi ~ 1.1 ± 0.3 Mq. Prom equation (2) 
and our adopted i = 45 — 75°, this implies the secondary 
mass is M2 ~ 0.6 ± 0.2 Mq, and the orbital separation 
a ~ 120 ±11 Rq. Therefore, the Roche lobe radii of the 
two stars are Rrls ~ 53±4 Rq, and Rrl,2 ^ 38±6 Rq, 
for the M giant and the secondary, respectively (Eggle- 
ton 1983). The Roche lobe filling factor of the M giant 
is then /rl = Ri/Rrls = 0.66 ± 0.17. 

If the deviation of the ASAS folded light curve and the 
ellipsoidal model with a giant alone is caused by a hot 
component due to a WD companion, its accretion disk 
and/or a tidal stream, then it is about 3% the luminosity 
of the giant in the V band. The luminosity of the hot 
component is then 2 — 8 L©, for an effective temperature 
in the range of 3800 - 10,000 K. For a M2 - 0.6 Mq 
WD with i?2 ^ 0.01 Rq, assuming all the gravitational 
potential energy is converted to thermal radiation, the 
corresponding accretion rate is 1 — 4 x 10^^ Mq yr^^. 

During the 1942-1950s outburst, the B magnitude of 
J0757 increased by ~ 1.5 mag at the peak. At a distance 
of D ^ 1 kpc. the absolute B magnitude of the outburst is 
MB,ouibursL peak ~ 0.7 mag. The peak bolometric lumi- 
nosity of the outburst, Loutburst peak dcpcuds on its SED. 
If the outburst has similar SED as a MOIII star, then 
Loutburst peak ~ 750 Lq. The minimum required bolo- 



metric luminosity would be when the radiation is peaked 
in the B band, i.e. Loutburst peak ^ 100 i© at T 8100 
K. Such a temperature is consistent with the observed 
hot component temperature of CH Cyg during its 1992 
outburst powered by accretion (Skopal et al. 1996). If 
the outburst is powered by Hydrogen shell-burning, the 
temperature during the outburst would be much hotter, 
i.e. Toutburst ~ 10^ K as observed in some symbiotic stars 
(Miirset et al. 1991). For example, Z And is found to 
have a temperature of ~ 9 x lO'' K during its 2000-2002 
outburst (Sokoloski et al. 2006). If the outburst in J0757 
has a temperature of 9 x 10^ K, most of the radiation is in 
the UV, and the corresponding bolometric luminosity of 
AB = 1.5 mag increase is Loutburst peak ~ 2.5 x 10^ 
assuming black body radiation. 

4. DISCUSSION 

The amplitude and timescale of the 1940s outburst, 
and its stability during quiescence phases (except the 
periodic ellipsoidal variation), is very unusual. The mass 
loss rate of a MOIII star is ~ lO'^ - 10"^ M© yr'^ 
(Cox 2000). With a Roche lobe filling factor ~ 0.5 - 
0.8 as discussed in section 3.2, a significant fraction of 
the wind could be transferred to its companion via wind 
Roche lobe overflow (Podsiadlowski & Mohamcd 2007). 
If the secondary is a 0.6 M© Main Sequence (MS) star 
with i? ~ 0.6 Rq, the only way to power the outburst 
is by gravitational energy release. The least required 
accretion rate to power the 1940s outburst of J0757 is 
then 3 x 10"*^ M© yr"! {Lp^ak ^ 100 L© if T - 8100 
K), which is 2 to 3 orders of magnitude higher compared 
with the typical mass loss rate of the M giant. Therefore, 
a compact object, almost certainly a WD (at ~ 0.6 M©, 
it cannot be a Neutron star or a black hole), is needed 
to power the outburst. 

If the outburst was powered entirely by accretion on to 
a 0.6 Mq WD, the least required accretion rate to power 
the 1940s outburst of J0757 is 5 x 10~* M© yr"\ assum- 
ing all the potential energy is converted into radiation 
{Lpeak ~ 100 L© if T ~ 8100 K). However, J0757 only 
showed one single outburst and was stable before and af- 
ter the outburst, which is different from other symbiotic 
systems with accretion-powered outbursts, such as PU 
Vul (before its nuclear runaway outburst in 1977), which 
was highly variable with multi-epoch outbursts (Liller 
& Liller 1979). Another symbiotic variable powered by 
accretion, CH Cyg, was stable prior to 1963 for several 
decades, but later showed a series of outbursts since the 
onset of activity in 1963 (Mikolajewski et al. 1990). 

The most probable explanation of the J0757 outburst, 
is then Hydrogen burning on the surface of the WD. 
The outburst profile of J0757, i.e. a relatively sharp 
rise and exponential-like decay, closely resembles that 
of Z And and CI Cyg, which are believed to be re- 
lated to the presence of unstable accretion disk around a 
white dwarf with hydrogen burning shell (Mikolajewska 
2003, et al. 2002). A "combination nova" outburst on 
Z And has been reported by Sokoloski et al. (2006), 
which is triggered by an accretion instability and re- 
sulted in increased thermonuclear burning. The outburst 
of J0757 may have started at the stable burning zone 
(M > Merit = 3.3 X 10-8 Mq yr-i for a 0.6 M© WD; 
Nomoto et al. 2007; Shen & Bildsten 2007) with an in- 
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TABLE 1 

ASTROPHYSICAL PARAMETERS FOR DASCH J0757. DETAILS ARE DISCUSSED IN SECTIONS 2Sz3. 



Parameter 

Observational parameters with highest confidence: 



Value 



Spectral type 
To (HJD) 
Period (days) 
7 (km s " ^ ) 
Ki (km s-i) 
Eccentricity 
uj (deg) 
Vrot sin i (km 



"Horn 

2455597.980 ± 0.054 
119.18 ±0.07 
-33.57 ± 0.04 
13.906 ± 0.025 
0.025 ± 0.01 
70.49 ± 0.09 
10 ±1 



Parameters derived from atmosphere fitting: 


^eff (K) 


3850 ± 50 


log g (cgs) 


1.0 ±0.5 


[Fe/H] 


-0.6 ±0.1 


Estimates based on models: 


i (deg) 


60 ± 15 


Ml (Mq) 


1.1 ± 0.3 


Ma (Mq) 


0.6 ±0.2 


Orbital separation (Rq) 


120 ± 11 


Ri (Re) 


35 ±9 


Li (Le) 


250 ± 130 


Distance (kpc) 


1.0 ±0.3 


hz (pc) 


400 ±110 


Mg in quiescence 


2.2 ±0.6 


Mb outburst 


0.7 ±0.6 


■^-outburst peak i^o) 


lO'^ — 10** (if no constraint on T/jot) 


100 (if Thot ~ 8100 K) 




2.5 X 10* (if That ~ 90,000 K) 


Lfiaf in quiescence (Lq) 


~2-8 


M in quiescence {Mq yr~^) 


~ 1 - 4 X 10"^ 



creased accretion rate due to an accretion disk instability 
(King et al. 2003), when the accumulated envelope mass 
on the WD exceeded the critical value (3 x 10~^ Mq; 
Townsley & Bildsten 2004). Then, consequently, hy- 
drogen was ignited and started quasi-steady burning to 
power the outburst. The luminosity of such a quasi- 
steady burning increases with higher WD mass, as well 
as higher accretion rate. For a 0.6 Mq WD, the theoreti- 
cal luminosity at McriL is L ^ lO^ ** Lq, and the effective 
temperature is 2.5 x 10'"' K; If M 4A'/crit. then the 
luminosity is ^ iO^ Lq, and the effective temperature 
is ~ 8 X 10'' K (Nomoto et al. 2007). Such a range of 
luminosity and corresponding temperature is consistent 
with the estimated peak luminosity of J0757, which is 
~ 10^ if T ~ 10^ K. 

The timescalc of the Hydrogen shell burning strongly 
depends on the mass of the WD, accretion rate, and 
whether a optically thick wind occurs (Kato & Hachisu 
1994; Prialnik & Kovetz 1995). Kato & Hachisu (2011) 
suggested that the companion star plays an important 
role in producing a wind, and thus affects the timescale 
of the outburst. As shown in Figure 6 in Kato & 
Hachisu (2011), in wide symbiotic binaries like PU Vul 
{Porb = 4900 days), the expanded WD envelope is not 
affected by the companion, resulting in a quiet evolution 
without strong winds in the outburst, and thus longer 
outburst (peak plateau timescale of ~ 10 yr); On the 
other hand, in close binaries such as V723 Cas and HR 
Del {Porb < 1 day) with similar WD masses as PU Vul 
0.6 Mq), the companion would have been deeply em- 
bedded in the nova envelope, resulted in strong winds 
and much shorter outburst (peak plateau timescale of 
- 200 days). We note that J0757, with Porb = 119.2 
days and a peak plateau lifetime of ~ 1 yr (as shown in 
the lower- left panel in Figure 1), lies interestingly in the 



region between the canonical symbiotic novae such as PU 
Vul, and novae in close binaries such as V723 Cas. 

The abscincc! of emission lines in the spectra of J0757 is 
peculiar, but not incompatible with the Hydrogen burn- 
ing powered outburst picture. Nova emission lines dis- 
appear quite rapidly after the outburst (Downes et al. 
2001), and it is not surprising that we do not see emis- 
sion lines 60 years after the outburst. The lack of emis- 
sion lines could be due to either or both of the following: 
The first possible reason is that there is no significant 
surrounding material. This may be either due to small 
amount of mass loss in the system even during the nova 
outburst, which is consistent with the recent results by 
Starrfield et al. (2012); or due to the fact that ejected 
shells disappear very fast, such as observed in V723 Cas 
(Krautter et al. 2002). The second possible reason is the 
lack of ionizing UV photons. A post-nova WD can cool 
off fast, for example, the white dwarf in DQ Her has been 
found to be below 10 Lq 40 years after the nova outburst 
(Stasifiska & Tylenda 1990). As discussed in sections 2 & 
3, the luminosity of the hot component in J0757 is only a 
few Lq, while the luminosity of the M giant is ~ 250 Lq, 
which dominates the spectrum. Similar cases have been 
found in a few symbiotic X-ray binaries, where the hot 
component is 3 orders of magnitude fainter than the 
red giant, and emission lines are absent in the optical 
spectra (Masetti et al. 2002, 2007). 

In summary, DASCH J0757 is an intriguing system 
with many unusual properties. It differs from other 
known classic or symbiotic novae, and is more likely a 
new type of nova-like variable located in the gap region 
between classic novae and symbiotic novae. Unlike most 
other symbiotic stars, it does not contain a luminous hot 
component or nebula emission lines, and its magnitude is 
stable after the 1940s outburst. Its current photometric 
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and spectroscopic profile is not different from a normal 
red giant binary. It would not be picked out without the 
detection of its long outburst in 1940s, made possible by 
the long timescales observed with DASCH. Future ob- 
servation in the UV and X-ray will help understand its 
mass transfer and accretion processes. 

Given the lack of other observations during the out- 
burst of J0757, it is hard to reach a definite conclusion 
whether it was powered by Hydrogen shell burning or 
accretion. If it was indeed Hydrogen burning with no 
significant mass loss, it could be a promising channel 
of accumulating mass on the WDs. Systems similar to 
DASCH J0757, but with more massive WDs (>- 1 Mq), 
could thus be promising type la SNe progenitors. They 
are not luminous supersoft X-ray sources in quiescence 
(see also Di Stefano 2010), and would not be picked out 
by the usual searches for symbiotic stars, but could be 
detected by the study of their long-term light curves. 
Note that J0757 is found from ~ 700 plates, while there 
are nearly half a million plates available to be scanned 
and analyzed in DASCH. With its 100 year light curve 
coverage, DASCH has a great potential to find more of 
such sources, which will shed light on the understanding 
of type la SNe progenitors, as well as the study of binary 
evolution and accretion physics. 
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Fig. 1. — Light curve of DASCH J075731. 1+201735. The upper panel shows the whole light curve. Black dots with error bars are the 
light curves from DASCH, and small green dots after Year 2000 are the light curves from ASAS. Since ASAS data are in the V band, while 
DASCH magnitudes are in the B band, we added 1.5 mag to the ASAS V mag in the plot to convert it to B, as described in the text. The 
lower left panel shows the light curve from 1940 to 1954 for the outburst event; Data before the red dashed line (May 1942) are used as 
before the outburst. The lower right panel shows the light curve around 1942 for the rising of the outburst. The blue dashed line marks 
the median magnitudes in JD 2430250 - 2430440 (12.22 mag; quiescence state) and JD 2430730 - 2430850 (10.68 mag; brightest state), 
and a naive linear increase in between; The time span of the linear increase is 88.5 days. 
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Fig. 2. — Folded light curve and period search of DASCH J0757. The upper panel shows the foled light curve with P = 119.18 days and 
To = 2455597.98. Black solid dots are DASCH data binned by phase, with error bars represent the standard deviation of the mean. Small 
green dots are ASAS data. Blue solid line is the best fit 3 harmonics Fourier model derived from the ASAS data. The lower panels show 
the reduced ^ ^ function of folding period. From the left to the right are the reduced using ASAS data, DASCH data excluding 
the outburst, DASCH data before the outburst (before June 1942, as indicated as red dashed line in the lower left panel in Figure 1), 
and DASCH data after the outburst (after 1960). The black dots are the reduced x~square values, and the red lines are a sixth order 
polynomial fit. In the two lower middle panels, i.e. DASCH data excluding the outburst, and DASCH data before the outburst, 1 — cr 
confidence levels by assuming distribution are marked as red horizontal lines. The more limited DASCH data after the outburst may 
suggest a ~ 0.02 days longer period, but this is not a significant deviation. The vertical blue dashed lines mark the best fit period of 119.18 
days. 
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Fig. 3. — FAST 300 grating spectrum of DASCH J0757. Spectral resolution is about 7 A. The upper panel shows the whole spectrum. 
The lower panels show two luminosity indicators, Ca I at 4226 A, and the CaH B-band at 6385 A with the blend near 6362 A consisting 
of Ti, Fe, and Cr on the left. 
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Fig. 4. — Radial velocity of DASCH J0757 as a function of orbital phase, using the ephemeris in Table 1. Blue solid dots are observation 
data using the TRES spectrograph at the FLWO 1.5in, and the smooth curve is the best-fitting model. Residuals from the fit are given in 
the lower panel. 



DASCH Peculiar Nova-like Outburst 



11 



A (A) 

10^ 10" 10^ 10^ 10^ 




10^' 10^' 10^' 10^' 10^' 

v(Hz) 



Fig. 5. — SED of DASCH J0757 using archival data. The black line represents a synthetic spectrum of a T^ff = 3750 K, and log g= 1.5 
giant (Lejeune et al. 1997). Note that SDSS riz bands are saturated; It was not detected by GALEX FUV and ROSAT. 
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Phase (P=1 19.14 d, Tpg^|=55, 593. 6794) 

Fig. 6. — Top: The phased and binned ASAS V-band Ught curve of J7057 and some representative eUipsoidal models. The red Hne 
is a model that includes irradiation and non synchronous rotation, and the blue line is a model that has an accretion disk with a hot 
spot. Bottom: The TRES radial velocity measurements phased on the same ephemeris (which uses the time of periastron passage as the 
zero-point) and the best-fitting model. 



